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O Abstract. This paper constructs a certain planar four-body problem which ex- 

hibits fast energy growth. The system considered is a quasi-periodic perturba- 
tion of the Restricted Planar Circular three-body Problem (RPC3BP). Gelfreich- 
Turaev's and de la Llave's mechanism is employed to obtain the fast energy 
growth. The diffusion is created by a heteroclinic cycle formed by two Lyapunov 
I I periodic orbits surrounding Li and L2 Lagrangian points and their heteroclinic 

intersections. Our model is the first known example in celestial mechanics of the 
a priori chaotic case of Arnold diffusion . 

a 

g 1. Introduction 

In this paper, we shall construct solutions of the Restricted Planar four-Body prob- 
^ lem (RP4BP) which exhibit long time instabilities, i.e. motions change substantially. 

CN The model we employ is that of a Sun- Jupiter-Planet (small)- Asteroid system. In 

this model the mass of the asteroid is assumed to be negligibly small (in fact zero). 
^ The mass of the planet, denoted mp, is strictly positive and < mp <C mj where 

^ mj denotes the mass of Jupiter. The Sun-Jupiter-Planet (S-J-P) system forms a 

T— I planar three-body problem (P3BP) which has quasi-periodic motion, and the ob- 

ject of study is to understand the motions of the massless asteroid in this system. 
Gelfreich-Turaev proposed the following mechanism for Arnold diffusion [2] . (It was 
_ ^ pointed out to me that R. de la Llave proposed the same mechanism in his unpub- 

^ lished paper [3] and mentioned it in his ICM 2006 talk.) For a Hamiltonian system, 

^ H{p, q, et), g G T"^, p G Mf^, consider the frozen system, H{p, q, u). Suppose for each 

u, and each energy surface of the energy interval of the frozen Hamiltonian 

H{p,q,iy), there are hyperbolic periodic orbits 71 and 72 with stable and unstable 
manifolds H^"'*(7i), ^^"'^(72) which make transversal heteroclinic intersections. In 
[2] adiabatic invariant theory and a special contracting mapping theorem are used 
to prove that, for a sufficiently small e, there exists ti > such that the Hamiltonian 
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Figure 1. Configuration of the four-body problem 

H{p,q,£t) lias a linearly fast diffusing orbit {p,q)(t), i.e. 

F(p(0),g(0),0) = h., H{p{t),q{t),et) = h+ 

for some t < ti/e. 

One important feature of the mechanism of pj is that their methods do not rely on 
how the Hamiltonian depends on et. The results hold true for periodic, quasiperiodic 
and other settings. The key to this paper is applying the mechanism to the 
RP4BP. In our case, the et dependence is quasi-periodic. 

Notice that, in Planar Restricted Circular three-Body Problem (RPC3BP), there are 
two normally hyperbolic periodic orbits 71, 72 surrounding the Li and L2 Lagrangian 
points respectively and their stable and unstable manifolds have heteroclinic inter- 
sections. This gives the "heteroclinic cycle" required in [2]. 

To obtain slow time-quasi-periodic perturbation we need to exploit the planet. We 
select S-J-P to have quasi-periodic orbit, so we shall have a quasi-periodic perturba- 
tion of the RPC3BP formed by S-J-A. The motions of the asteroid can be described 
as solutions to Hamilton's equation with a Hamiltonian Ha that is a quasi-periodic 
perturbation of the RCP3BP. (see equation Q): 

HAiLA,iA,GA,9A,t) = RPC3BP+6fiLA,£A,GA,9A,et), {LaJa,Ga,9a) G T*T\ 

Here we have two independent small parameters, e is the frequency of the S-J-P 
system and 6 is determined by the mass of the planet mp. The variables we are 
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using here are called Delaunay coordinates. See Section 3 and Appendix [X] for more 
details. The main theorem proved in this paper is: 

Theorem 1. In the Sun-Jupiter-Planet- Asteroid system, for small e and 6 > 0, 
there exists a diffusion orbit {La, ^a, Ga, 9A)it) with linearly fast speed energy growth, 
i.e. There is an energy interval [h^,h^], such that the energy of the asteroid has 
growth: 

HAi{LA,iA,GA,gA)iO),0) = /i-, HAiiLA,iA,GA,gA){t),et) = h+ 
for some t < ti/6e, for some ti, where /i_ are independent of 6 and e. 

The problem of the Arnold diffusion is a long story concerning the instability of 
generic Hamiltonian systems. Here we do not try to mention seas of literatures 
about Arnold diffusion, but point out the results relevant to our work. Moreover, 
even though the problem of Arnold diffusion has been studied for half a century. 
There are scarce concrete examples, esp. in celestial mechanics. As far as the 
author's knowledge, the only known examples are |14lll5t [T6 t ll7]. Their mechanisms 
are close to the spirit of Arnold's example, called a priori unstable case. 

But our model has new features. The study of energy growth is a simplified version 
of Arnold diffusion by Mather, so it is also called the Mather problem [T]. The 
mechanism of diffusion is called "a priori chaotic" by [1], since the reference system 
has some conserved quantities, but there are orbits which are hyperbolic and with 
transverse heteroclinic intersections in the manifolds corresponding to the conserved 
quantities. The systems are not close to integrable, the Nekhoroshev upper bounds 
for the time of diffusion does not apply. Our model is the first known model of the 
a priori chaotic case in celestial mechanics. It turns out we obtain linear energy 
growth rate in our system. 

The paper is organized as follows. In Section 2, we give a brief introduction to the 
Gelfreich-Turaev's mechanism. 

In Section 3, we give the construction of the configuration of the four-body problem. 
We first find a family of quasi-periodic orbits for full S-J-P three-body problem. 
Here we use Poincare Continuation method. Then we fix the quasi-periodic orbit 
and write down the Hamiltonian governing the motion of the asteroid as a quasi- 
periodic perturbation of the RPC3BP (S-J-A). 

In Section 4, the RPC3BP (S-J-A) is studied. There are two normally hyperbolic 
periodic orbits 7i(/i), 72(^) around the Li and L2 Lagrangian points respectively, on 
each energy level h for an energy interval h G These periodic orbits 

71 (/i) and 72 (^) are the key ingredients for the application of Gelfreich-Turaev's 
mechanism. It is also shown by computer assisted proof in [18] that the stable 
and unstable manifolds of the normally hyperbolic periodic orbits have transversal 
heteroclinic intersection. 
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In Section 5, the heteroclinic cycle of the RPC3BP (S-J-A) is transplanted to the 
RP4BP. This is done using the hyperbolic theory. 

In the last section, Gelfreich and Turaev's mechanism is applied to the RP4BP. 
With the heteroclinic cycle established in Section 3,4,5, the uniformity assumptions 
[UAl] [UA2] and the nondegeneracy condition required in [2] are verified. So the 
main theorem follows. 

2. A BRIEF INTRODUCTION TO GeLFREICH-TuRAEV'S MECHANISM 

Consider a Hamiltonian system H = H{p,q,et), {p,q) G M^" with e small. It is 
routine to consider the frozen system in adiabatic invariant theory H = H{p, q, u), 
where = et is treated as a constant parameter. It is required that the frozen 
system has a chaotic behavior, namely there exists uniformly-hyperbolic, compact, 
transitive, invariant set A^u in every energy interval H £ [h-,h^] for all u. In 
every given energy level, the set A^,^ is in the closure of a set of hyperbolic periodic 
orbits each of which has an orbit of a transverse heteroclinic connection to any of 
the others. This means that orbits of the frozen system may stay close to any of 
the periodic orbits for an arbitrary number of periods, then come close to another 
periodic orbit and stay there, and so on. Now we take two periodic families 71 and 
72 of the frozen system. It is shown that under some natural conditions one can 
arrange jumps between 71 and 72 in such a way that the energy keeps growing. It 
is proved that 



where Ti is the period of the periodic orbit 7j. Assume that the differential equation 



has a solution hs{i^) for 6' sufficiently small to suppress j3 ( where (3 is defined 
in equation (46) 0/ ^). Assume the uniformity assumptions [UAl] and [UA2] 
hold true. Then for all sufficiently small e the Hamiltonian system has a solution 
{p{t),q{t)) such that 



We shall show hs grows linearly. Note that the uniformity assumptions , [{7^42] 

are automatically fulfilled for any compact set of h and v, which is exactly what we 
consider. So we do not cite the lengthy formulation of [C/Al], [[/A2]. 



Theorem 2 (Theorem 2 and 3 of [2j). Let 




dh 
du 



max{t)i(/i, I/), U2(/i, v)} — S' f3{h, v) 



H{p{0),q{0),0) = hs{0), H{p{t),q{t),et) = hs{et). 
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3. The configuration of the four-body problem 

In this section, we first establish the quasi-periodic motion of the Sun- Jupiter-Planet 
system, then write down the Hamiltonian governing the motion of the asteroid. 

3.1. Selection of quasi-periodic orbits of the P3BP formed by Sun- Jupiter- 
Planet. We need to have a family of quasi-periodic orbits for P3BP (S-J-P) with 
the following properties: 

(1) Orbits of Sun and Jupiter are nearly circular. 

(2) Orbits of the planet are nearly elliptic, i.e. osculating eccentricity is nearly 
constant. 

The reason that we need a "family" of quasi-periodic orbits is only technical. We 
need this family to check the nondegeneracy condition in the last section. 

We first show the existence of periodic orbits in the S-J-P system when nip = 
0, rrij 7^ 0, and Sun-Jupiter has circular orbit, i.e. the RPC3BP. Then we continue 
it to the case of mp > 0. 

Before the proof, let us collect the notations in the following definition. Different 
coordinates will get involved for the convenience of proofs, such as the Cartesian co- 
ordinates and Delaunay coordinates- • • . Please go to Appendix |X] for the derivations 
and physical meanings of them. 

Definition 1. (1) In the Cartesian coordinates, we use {x,x,y,y) {or {q,p)). 
{x,y){ or q) is the position and {x,y){or p) is the velocity. 
(2) In Delaunay coordinates, we use the variables {L,£,G,g). iP' is the semima- 
jor of the Keplerian ellipse, i the mean anomaly, G the angular momentum 
and g the argument of the perihelion. 

Definition 2. (1) We use the subscript to indicate the corresponding quantity 
of a certain body. For example, Gp, the angular momentum of the planet, 
rp, the mutual distance from the planet to the origin. 
(2) rpj,rps,rAj-,rAS^fAP denote the mutual distances between the planet and 
the Jupiter, the planet and the Sun, the asteroid and the Jupiter, the asteroid 
and the Sun, the asteroid and the planet, respectively. 

In the following theorem, we prove the existence of a family of periodic orbits for 
the RPC3BP. 

Theorem 3. The P3BP{S-J-P) has a family of quasiperiodic orbits. This family 

\ 1 + al 

can be parametrized by the eccentricity of the planet Cp G C, 1 where C is 

ap 
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such that the length of the interval 



C, 1 - 



is 0(1) as mp,mj — )• 0. When 



ap 
1 + a 

ep = 1 , the closest distance between the planet and Jupiter is a + o{mp) as 

ap 

mp — )• 0. Each of these quasiperiodic orbits covers a torus T^. 
Proof. The proof is split into 3 parts. 

1) Set mp = to find one periodic orbit in the RPC3BP (S-J-P). 

2) Find a family of periodic orbits in the RPC3BP (S-J-P). 

3) Continue these orbits to full 3 body problem (mp 7^ 0). 
STEP 1, find one periodic orbit in the RPC3BP (S-J-P) 

In this step, mp = 0. The motion of the planet is a RPC3BP (S-J-P). The Hamil- 
tonian for it in Delaunay coordinates and rotating coordinates is given by (c.f. [7J, 

my- 

Hp = L - Gp + Ai/, where A// = — - ^ - L^. 

2L|, rp rpj rps 

Here we write the perturbation using the polar coordinates for simplicity. It should 
be converted into the Delaunay variables. Rotating coordinates is the noninertia 
coordinates we choose to fix the Sun and Jupiter on the x-axis. 

The plan is to consider the double Kepler problem (2BP (S-P) -|- 2BP (S-J)) by 
setting mj = 0, then the Hamiltonian in rotating coordinates becomes 

1 

iUi 



Hp = -77T2" ~ ^P- 



Consider the resonance relation gp = Nip, where N = p/q, p,q £ Z is a very large 
rational number. This implies, after p period of the planet, Jupiter experiences q 
periods. The whole S-J-P system comes back to the original configuration. In the 
rotating coordinates, this is a periodic orbit for the double Kepler problem. The 
number N is large so that the period is long. We want to find a periodic solution 
of the RPC3BP (mj > 0) in a neighbourhood of this resonance. If is a large 
number, then the semimajor of the ellipse of the planet should be very large. It is 
proven by Arenstorf and Barrar that the periodic orbits that are symmetric along 
the X-axis are locally isolated and therefore can be continued to the RPC3BP for 
mj > 0. 

Lemma 1 (c.f. [9]). Suppose when mj = 0, yp(0) = and xp(0) = 0, i.e. at time 
t = 0, the planet crosses the x-axis perpendicularly, and for a later time T/2, the 
planet again crosses the x-axis perpendicularly, i.e. yp(T/2) = 0, xp{T/2) = 0. 
Then except at most finitely many eccentricity ep, there exists a periodic orbit of 
the RPC3BP(S-J-P) for each ep and for mj > sufficiently small, whose period is 
slightly different from T, with error o{mj). 
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We use the Keplerian relation 



1 



that relates the semimajor a to the period 



r2 47r2 

T. Moreover, we have a = relating the Delaunay variable L to the semimajor. 
Since we consider fixed N, we fix the period T and hence the semimajor ap and Lp. 

STEP 2, find a family of periodic orbits for the RPC3BP (S-J-P) 

To find a family periodic orbits for the RPC3BP (S-J-P), we notice in the Kepler 
problem, the period of the ellipse depends only on the semimajor. Varying the 
eccentricity, we obtain a family periodic orbits parametrized by ep in the double 
Kepler problem without violating the resonance relation gp = Nip. For each one of 
them, we can use Lemma [T] to continue it to positive mj. Now we want the closest 
distance between the planet and Jupiter to be a > 0. Let us see how large should 
ep be. The orbit of Jupiter is a circle of radius 1 centered at the Sun and that 
of the planet is an ellipse focused at the Sun. From geometric consideration, we 
immediately obtain 

ap — apep > 1 + a. 

I + a n \ 1 + a 

This gives us ep < 1 . Now we apply Lemma ItoepG 0,1 with 

ap 



possibly finitely many exceptional Cp's. We choose some C G 



0,1 



ap 
1 + a 
ap 



such 



that the interval 



C,l 



1 + a 
ap 



uniform mj for all the ep because of the compactness of 



contains no exceptional ep's. Then we can find a 

1 + a" 



C,l 



ap 



STEP 3, continue the family of periodic orbits of the RPC3BP (S-J-P) to the P3BP 
(S-J-P) 

Now we use Poincare's continuation method to continue the periodic orbits in the 
RPC3BP to quasi-periodic orbits in the full 3 body problem. The result was first 
proved in p!0] , and the proof is reformulated in a concise form in [8] . 

Lemma 2. (Theorem 9.6.1 in [8j) 

Any elementary periodic solution of the planar restricted 3 body problem whose period 
T is not an integer can be continued into the full 3 body problem with one small mass. 



Here, "elementary" means, the periodic orbit is "isolated" on the energy surface. 
Our periodic orbit found in Step 1 has this property. 

In order to apply the theorem in our case, we only need to ensure that the period 
T is not an integer. In fact, We get a periodic orbit of the RPC3BP (S-J-P) from 
a Keplerian elliptic orbit. In order to avoid the danger that T G Z, we allow in the 
resonance relation gp = Nip that is a fraction. After continuing a periodic orbit 
of the RPC3BP to the full 3BP(S-J-P), the orbit is no longer periodic in general 
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but is quasi-periodic. The dimensional consideration shows that this quasiperiodic 
orbit covers a torus T'^. □ 



3.2. The motion of the asteroid. In this part, we consider the motion of the 
asteroid under the force from the Sun, Jupiter, and the planet. The following 
theorem establishes that the Hamiltonian describing the motion of the asteroid in 
the 4BP(S-J-P-A) is quasi-periodic perturbation of the RCP3BP(S-J-A) with 2.5 
degrees of freedom. This is the Hamiltonian for which we prove the existence of 
diffusing orbits for this Hamiltonian. 

Theorem 4. // the motion of S-J-P is chosen to be quasi-periodic, the Hamiltonian 
of the motion of the asteroid can be written into the form: 

1 



(1) HA{LA,iA,GA,gA,£t) 



2L\ 



Ga + AH + 5fiiA, LA,gA,GA, et), 



where {LA,iA,GA,gA) G T*T^, the 6 = o{mp) as mp 0, f is quasi-periodic w.r.t. 
t, and e = 1/N for some large N. i.e. Ha =RPC3BP-h quasi-periodic perturbation. 



Proof. The Hamiltonian of the motion of the asteroid can be written in the complete 
form: 



HA{xA-,XA,et) = ]^x\ 



1-mj 



TAJ 



mp 



{xa,xa) G 



TAS rAj rAP 

The first term is the kinetic energy and the last 3 terms are the potential energy 
from the Sun, Jupiter and the planet respectively. 

We want it to be close to the Hamiltonian of the RPC3BP (S-J-A). The problem 
is the orbit of the Jupiter is no longer circular, because the mass of the planet is 
positive. We first choose the mass center of the Sun, Jupiter and the planet as the 
origin, then write the Hamiltonian in rotating coordinates. We want that, in the 
rotating coordinates. Sun and Jupiter lie on a line parallel to the x-axis. As a result 
this rotating coordinates is not uniform. We use the following transformation. 

\ qA -= exp(£jir) XA 
I PA ■■= exp{ejK) XA 



where K= 



1 

-1 



and £j is the mean anomaly of the Jupiter after taking into 
consideration the perturbation by the planet, see Appendix [A| equation (12). 
exp{ijK) is an orthogonal matrix. So it is easy to check the following identity. 
jK n 1 r n 1 r p£jK 











Id 
-Id 









ijK 



1 ^ 


" Id' 




-Id 
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This shows that the change of coordinates is symplectic. The coordinate system 
(qAtPa) is also Cartesian. 

ij 7^ constant, so the rotating angular velocity of the coordinate frame is not 
uniform. Instead of —q^KpA, the term 

(OH \ 
1 + "ip^^j QaKpa 

would appear in the new Hamiltonian as the Coriolis term (c.f. the 6th chapter of 

|8]), where mpHi is defined to be the term — — - — ^ — J—, i = J, 2 = P in 

mo IQ1-Q2I 



equation (11) in Appendix B I Now we get the Hamiltonian: 



Ha = 7,Pa - ijQAKpA 

2 rAS rAj TAP 

Sun and Jupiter are very close to the points (— mj,0), (1 — mj,0). This deviation 
was explained and estimated in Lemma 3.1.3. Now plug in the terms: 

1 — mj nij 1 — mj mj 
n r2 ri r2 

(ri and r2 are defined to be the distance from the asteroid to {—mj, 0), (1 — mj, 0) 
respectively. The difference of ri to vas and r2 to r^j is o(mp).) 



So the Hamiltonian becomes: 



Ha 



1 2 Tr^ '^-mj mj 
\Pa - QaKpa 



+ 



dHi Tj^ 
mp-^QAKpA 



2 ri r2 

+ (l-mj)(-- — ) — ) - — 

In this expression, the first bracket is the RPC3BP in a uniform rotating coor- 
dinates. The second bracket is of order o(mp) since it is zero when mp is zero 
and the Hamiltonian system depends on mp analytically. When we plug in the 
information of the motion of Sun-Jupiter-Planet {£j, Lj, gj,Gj, ip, Lp, gp,Gp) as 
a quasi-periodic vector-valued function of t. The second bracket (denoted by df) 
depends on time t quasi-periodically. 

The Hamiltonian Ha is written in terms of the Cartesian coordinates {qA,PA)- We 
can convert the Cartesian coordinates into Delaunay coordinates to obtain ([T]), where 
the first 3 terms are the Hamiltonian of the RPC3BP and 
(2) 

dHi J, \( ^ ^ \ , ( ™p 



(5/(£A,iA,ffA,GA,ei) := -mp— — g^i\:j5A + (l-m,7) — 

oLj V ri TAs \ r2 taj J tap 



Here the 5 = o{mp) as mp — )• 0. We use e = 1/A^ to denote the frequency of the 
periodic orbits that appear in the proof of Theorem [3] when mp = mj = 0. □ 
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Remark 1. For the perturbation we have 

(1) the term q^KpA blows up when Hi tends to infinity, i.e. Jupiter collides 



1 1 \ /I 1 \ 



dLj 
with the planet. 

(2) the term (1—mj) ( ) +mj ( ) blows up when each of the 

\ri rAs J \r2 taj J 

denominator tends to zero. i.e. the asteroid collides with the Sun or Jupiter. 

(3) the term — — blows up when rAP tends to zero, i.e. the planet collides with 

rAP 
the asteroid. 



4. Existence of heteroclinic cycle for the RPC3BP (S-J-A) 



In this section we neglect the influence of the planet in the 4BP(S-J-P-A) as defined 
in above in Section [3j The goal is prove the existence of two normally hyperbolic 
periodic orbits and the existence of transverse heteroclinic connections between them 
for the RCP3BP(S-J-A). In Section[5| it will be shown that these phenomena persist 
to the 4BP(S-J-P-A) when we let mp > 0. 



4.1. The existence of two hyperbolic periodic orbits in the RPC3BP. We 

describe the periodic motions 71 (resp. 72) of the asteroid near the Lagrangian point 
Li (resp. L2). Here, it is more convenient for us to study the motion in (x,y,x,y) 
Cartesian coordinates. In this coordinates (rotating coordinates), the equations of 
motion are (c.f.[7J): 

. n x'^ +y^ 1-mj mj mj{l - mj) 
where \ l = 1 ! 1 , 




I ri r2 I 

and ri = ■\/(x + vtljY + y'^ , r2 = sj (x ^ mj - 1)2 + i/2_ 
It is easy to check that these equations have an integral — the Jacobi integral. 

(3) J{x,y,x.y) := -{x'^ + f) + 2n = -2h, 

where h is the energy of the RPC3BP. The RPC3BP has 5 equilibria called La- 
grangian points (see Figure 2). We want to pay attention to the collinear equilibria 
Li,L2 lying on the x-axis. The Li and L2 Lagrangian points are the two positive 
critical points of the potential O restricted on the j;-axis (c.f. Chapter 2.5 of [7J): 

x'^ 1 — mj mj 

(4) + 1 — , ^ + 



2 Ix + mjl |x + mj — 1| 
We immediately have the following lemma 

1 /3 

Lemma 3. The distance from Jupiter to the Lagrangian points Li and L2 is 0{m j ) 
as mj — >• 0. 
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Figure 2. Five Lagrangian equilibria of the RPC3BP 



Proof. We consider L2 only, i.e. x > I — mj. The Li case is similar. Then Q 
becomes 

1 — mj mj 



+ + 



2 X + mj X + mj — 1 
Taking derivatiave and setting the derivative to be 0, we get an equation 

1 — mj rrij 



(x + mj)^ {x — 1 + rrij)^ 



0. 



As a first step approximation, we suppose x = 1 — mj + Cnij for some constant C 
and a and plug it into the equation to get 

1 - mj + Cm}(l - mj)(l - 2Cm'}) - Cm]-'' = 0. 

This equation holds only if a = 1/3. □ 

Since Li and L2 are critical points of the potential Q, we can linearize the Hamil- 
tonian system in a neighbourhood of the two points. 

Lemma 4. ([3]) The linearized systems in a neighbourhood of Li and L2 have eigen- 
values Ai^2 = iA, A3^4 = ibm, the corresponding eigenvectors, ui,U2,wi,W2, and the 
general solution: 



(5) u{t) = aiuie^* + a2W2e"^* + 2Re{^e'''^wi). 
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From this Lemma, we can see, for the hnearized system, the two complex conjugate 
purely imaginary eigenvalues give rise to periodic orbit, and the two real eigenvalues 
give the stable and unstable directions. 

It is proved in [3] that: 

Theorem 5. ([3]) The periodic orbits that correspond to ai = 02 = 0, project into 
the X — y-plane as an ellipse with major axis in the direction of y-axis, and minor 
axis in the direction of the x-axis . The direction on the orbit is counterclockwise, 
{see Figure 3) 

The proof of the theorem is an application of the Lyapunov center theorem (c.f. [3] 
and [8], Theorem 9.2.1). 

4.2. The existence of heteroclinic intersection in the RPC3BP. To prove 
the presence of transversal intersections for Vl^*(7i(J)) and M^"(72(^)), and vice 
versa, we cite the result from [18j. In [18j, the authors study the Hill problem that 
is a simplification of the RPC3BP in the limit mj — )• 0. Since nij is small, we 
consider the Hill limit mj — t- with proper translation and rescaling. Namely, we 
consider the following (c.f. [18J). 

X = Xruj^ + mj — 1, X = XruJ^ , y = Yrn^J^ , y = Yruj^. 

The rescaling is natural in view of Lemma [3} Then we express J in terms of 
{X, Y,X,Y), to obtain the Jacobi constant in the Hill's limit 

JniX, Y, X, Y) = -X^ -Y^+ ^^^i/^ + + 0{w}f'), 
where Jh = mj^^^{J — 3(1 — mj)) and rn'J^ Jh —s- J — 3 as mj — 0. 

It is proven in p!8] that 

Lemma 5 ([E])- For the Hill problem, for -^Jh\~'^^'^ > 1/18 the invariant mani- 
folds of the Lyapunov periodic orbits give rise to homoclinic and heteroclinic orbits, 
connecting a vicinity of Li to itself and to one of L2 and vice versa. These results 
extend to the RPC3BP for mj sufficiently small. 

This lemma gives us the heteroclinic intersections between the invariant manifolds 
of the Lyapunov orbits for RPC3BP(S-J-A) for sufficiently small mj. 

5. Existence of heteroclinic cycle for four-body problem 

In the previous sections, we have shown the existence of the heteroclinic cycle for 
the RPC3BP. Now we want to establish the structure for our four-body problem. 
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According to the Hamiltonian Ha ([T|, our four-body problem is a perturbation to 
the RPC3BP (S-J-A). The RPC3BP (S-J-A) and P4BP (S-J-P-A) have different 
degrees of freedom and phase spaces. So first we use the classical hyperbolicity 
theory to show the persistence of the two normally hyperbolic invariant cylinders 
and persistence of heteroclinic intersections under small perturbation. Second, we 
show that after freezing the P4BP the perturbation does not destroy the periodic 
orbits within the two cylinders. 

5.1. Translation of Section [4] into the language of classical hyperbolicity 
theory. Recall that the RPC3BP has 2 degrees of freedom and its phase space 
is r*T^. Denote the flow of the RCP3BP by ^>f. Fixing the energy restricts the 
dynamics to a 3 dimensional energy surface denoted by M{h), h G (In the 

following sections, we use h instead of the Jacobi constant J to stand for the energy 
to be consistent with the main theorem, J = —2h.) On every energy surface M{h), 
there are 2 hyperbolic periodic orbits 7i(/i), 72(/i) (see Theorem [s]). 

Definition 3. Define the cylinder Ci = \JfiJi{h), i = 1,2. These cylinders are 
two dimensional, and each is diffeomorphism to x T-*^. Additionally their 

stable and unstable manifolds, denoted W^{Ci) and W'^{Ci) respectively are both 
three dimensional. The heteroclinic intersections of these manifolds 

ri2 = [jri2{h) c w'iCi) n w^{C2), = (jT2i{h) c w^"(Ci) n w'{C2) 

h h 
are two dimensional objects. 

The following lemma is a translation of Section [4] into the language of classical 
hyperbolicity theory. 

Lemma 6. ([IJ) For some constant K, A > {calculated at the end of Section 4-1 
), and for all x G ^i{h), we have the decomposition of the tangent space into stable, 
unstable, and central subspaces. 

T^M{h) = E'^eE^eT,-fi{h) 

with 

\\D'^t{x)\E^J\ <Ke-^\ fort>0, 
\\D<l>t{x)\E^\\ <Ke^\ fort<0, 

ll^^t(^)lT.7,wll for tern. 

The stable and unstable manifolds to Ci: W^^'^\Ci), are three dimensional manifolds 
diffeomorphic to [h- , x T-"^ x M, and their heteroclinic intersections T12 and 
are both diffeomorphic to x M. 



IINXIN XUE DEPARTMENT OF MATHEMATICS, UNIVERSITY OF MARYLAND, COLLEGE PARK, MD, 20742 



5.2. Add the perturbation to the RPC3BP (S-J-A). When the perturbation 
is added, we start to consider the R4BP. The Hamiltonian is ([T]): 

Ha = -^-Ga + AH + 6f{£A, La, qa, Ga, et) := Hq + 6f, 

where Hq := 2 — Ga + Aif is the RPC3BP (S-J-A). The system as stated is 

non-autonomous. To apply the mechanism of [2j, the frozen the system with et = 
u, u gT'^ is considered. Note the frozen system is autonomous and its Hamihonian 
can be written in the form: 

HA{iA,LA,gA,GA,i^) = Ho{£A,LA,gA,GA) + 5 f {£a, La, gA, Ga, 1^) ■ 

Now we are working with an autonomous system with 2 degrees of freedom. For 
each fixed € T^, using the hyperbohcity theory we have the foUowing lemma. It 
shows that when 5^0 the cylinders along with their stable and unstable manifolds 
and heteroclinic intersections persist. 

Lemma 7. Assume the Hamiltonian has the form Ha = Ho + Sf^v), u gT'^ fixed, 
Ha G C"',2 < r < 00, Then there exists a 6^, such that for \5\ < 6^, the perturbed 
cylinder Ci^s is hyperbolic, G^~^ diffeomorphic to Ci, locally invariant and is 6-close 
to the unperturbed cylinder Gi. Its (un)stable manifolds W^^'^\Cs,i) is also 5-close to 
the unperturbed W'^^'^\Gi) in the C"""^ sense. (i=l,2). Moreover, the compactness 
o/T^(9 v) gives us a uniform 5* > such that the above statement holds for all 
1/ G T2 and \5\ < 6*. 



Proof. Direct application of classical hyperbohcity theory, c.f. pj, theorem 4.2 and 
theorem A. 14, A. 12. □ 

Remark 2. (1) Since we have the transversal heteroclinic intersections T12 C 
W'{Gi) n T^"(C2) and T21 C WiGi) n W%C2), there exists locally unique 
new heteroclinic intersections: Ts^i2 C ^^^(C^^i) H W^{Gs^2) o.'ri'd Ti^2i C 
W''{G5^i)r\W'{G5^2)- r5,i2 {Ts,2i) is 5-close toTi2 {T21) in the C^-^ sense, 
and that Ts^i2 (r5^2i) can be parametrized by a C""~^ function on T12 (J^2i) 
to the extended phase space. 
(2) The length of the cylinders, i.e. the energy interval may become 

shorter after the perturbation. Their difference would be of order 5, which is 
very small compared with — so we neglect it in the following discussion 
for simplicity of notations. 

We still need to show the Lyapunov orbits cannot be broken by the perturbation 
5f{v). Then we do not have big gaps within the cylinder. 

Lemma 8. Within each perturbed cylinder, the periodic orbits do not break if we 
add the perturbation 6f to RPC3BP. 
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Proof. Because our system is frozen, for fixed z^, it is a system of 2 degrees of 
freedom. As a matter of fact, the persistence of these perturbed periodic orbits 
can be estabUshed from the Lyapunov center theorem ([8], Theorem 9.2.1) in the 
same way as the existence of Lyapunov orbits. The Lyapunov orbits estabhshed 
in Theorem [5] is a perturbative result from the hnearized system Lemma |4| If 
we treat the 5f{u) and the nonhnear part of the RPC3BP expanded at Li or L2 
Lagrangian points as a whole, to perturb the linearized system, we still have a family 
of Lyapunov periodic orbits due to the same Lyapunov center theorem used by H] 
and [8j provided 5 is small. Since we have S T^, the compactness of gives a 
uniform 5 which works for all z^'s. □ 



6. Apply Gelfreich-Turaev's Mechanism to the four-body Problem 



6.1. Verification of the uniformity assumptions in [2]. Notice the P4BP under 
consideration satisfies the Gelfreich-Turaev's mechanism exactly. Indeed we have 
shown the existence of two normally hyperbolic periodic orbits and their transversal 
heteroclinic intersections on every energy surface of the frozen system in the energy 
interval In order to apply Gelfreich-Turaev's result, we need to check the 

uniformity assumption [UAl], [UA2] in [2]. We do not want to cite them here. 
[UAl] is the hyperbolicity requirement, which is given by Conley's result in Section 
4.1, while [UA2] is trivially satisfied, since we only consider finite energy interval 
which has compactness. Now we have proved the Theorem [l] modulo checking 
the nondegeneracy condition. The time span ti/5e can be seen from the following 
equations ([6]), ([s]). 



6.2. Verification of nondegeneracy. We break the proof into several steps. First 
we write the nondegeneracy condition into a form that we are able to check. Then 
we show the resulting form is analytically dependent on the variables. Then we use 
the property of analytic functions to show the nondegeneracy holds. 



6.2.1. Write integrals responsible for non- degeneracy. In the Hamiltonian (jl]), we 
set et = V. According to Theorem [2| we have 

(6) - = max [^^^,^^^\ - ^ m 
where 

1 f'^' 

(7) fi = — f\^^{lA,LA,gA,GA,iy)dt, 



and Ti (i = 1, 2) is the period of 7^. The 5' here is the 5 in [2|, which is just a small 
number, much smaller than 5 in this paper. So the 5' (5 term is neglible. 
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Now, integrate both sides. 

(8) 

h{h)-h{0) = l<5(/i(ti) + /2(ti)-/i(0)-/2(0))+<5^ ' 



ifi - h) 



du-25' / ^du. 



So in order to get linear energy growth, our non-degeneracy condition is (Theorem 
4 in [2]): 

(/1M-/2M) 



1 

hm inf — / 



du 



dv > 0. 



In our case, the perturbation is quasi-periodic. So it is sufficient to satisfy: 
(9) h{v) - Uu) ^ const. 

This is also the equation (69) in |2]. 

6.2.2. These integrals depend analytically on parameters and variables. We show the 
following: 

Lemma 9. The integrals ([T]) appearing in Q depend analytically on all the pa- 
rameters. 



Proof. Originally, away from any collision, the perturbation obtained from Sec- 
tion 3.2 the perturbation 5f (equation ([2])) is an analytic function in terms of the 



positions and momenta of Sun, Jupiter and the planet because of the analyticity of 
the Newton potential. 

For the S-J-P system, we can transform the variables into Delaunay coordinates: 

(^j, Lj, gj.Gj, £p, Lp, gp, Gp). 

The transformation is analytic, so / is analytically dependent on the Delaunay co- 
ordinates. 

Moreover, once we have chosen a quasi-periodic solution of the S-J-P system, these 
Delaunay coordinates would be a quasi-periodic function of time t. These quasi- 
periodic functions are also analytic because they are the solutions of the Hamiltonian 
equations, which are analytic ODEs (here we use Cauchy's theorem). 

So we get the Delaunay coordinates of the S-J-P system are analytically dependent 
on the time t. □ 



6.2.3. Find terms which should dominate these integrals and become unbounded. 
Lemma 10. Consider fi and f2 defined in we have 

/i(i^) - /2(i^) / const, 
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for some e p satisfying 
ofmj,mp. 



ep- 1 



1 

I2 



> c where the constant c > is independent 



Proof. We have shown that /i(z^) is analytically dependent on v and ep. Now 
consider the Fourier expansion 



/.(z.) = ^ Fkiep)e''\ V G T^. 

To prove the theorem, we only need to show Fk{ep),k G \ {0} does not vanish si- 
multaneously. We know F^s dependence on ep is analytic since ep = — [Gp/LpY 
Let us first fix some ep and vary v. We want to show the function fi{v) has a sin- 
gularity for some while at some other point v, it is bounded. This implies, for 
this fixed ep, fi{y) is not a constant function of v. So Fk{ep), k G I? \ {0} does not 
vanish simultaneously. Each Fk{ep) is analytic, so the zeros are isolated. We only 
need to avoid finite possible values of ep, e.g. the zeros of -^{o,i)(ep) to make fi{y) 
have nonvanishing dependence of u. 



We have found a family of quasi-periodic orbits of S- J-P which can be parametrized 

1 + a] n 

C, 1 :py — in Theorem 3[ If we choose a = dist[L2, (1 — mj, 0)) = 



by ep e 



7^2 



0{rnj^) according to Lemma 



then we expect when ep 



1 



1 + a 



the planet 



passes close to L2 and is at a distance of 0{mj'') from Jupiter. We are going to 
show that for such ep, the expression /2(i^) — /i(z^) blows up. Therefore using the 
properties of analytic function, we get /2(i^) — fii^^) is not a constant for any choice 
of ep. In the end, we choose ep = C. This configuration of P4BP (S-J-P-A) will be 
away from collisions and give us the desired linear energy growth. 

Now let us study the function 5f and its averages 5fi. We first consider the term 
■ 1 1 \ /I 1 ' 



(1 — mj) ( ) -|- mj ( ) in (|2|). We claim that the contribution of 

Vn rAs J _\r2 TAJ J 1_ 
this term to the averages 6fi is o{mp) as mp —t- 0. Indeed, we have the estimates 

1 1 \ Ar^-rAR\ . , dist{qs, {-mj,0)) 



(1-mj) 



rAS 



(1-mj) 



rirAS 



nrAS 



o{mp) 



as mp —>• 0, since the distance between the sun and the point {—mj, 0) is 

1/3 

dist{qs, {—mj, 0)) = o{mp) as mp — and ri,rAs = 1 — 0{m j ) as mj —)■ 0. We 
also have 



mj 







_ 'mj\r2 - rAJ 




rAJ J 


r2rAJ 



< 



mjdist{qj,{l-mj,0)) 2/3. 
= o(rnp)0{m J ) 



r2rAJ 
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as mj,mp — )• 0, since the distance between the Jupiter and the point (1 — mj, 0) is 

1/3 

dist{qj, {—mj, 0)) = o{mp) as mp — )• and r2, r^j = 0{m J ) as mj — )• according 
to Lemma [3l 



Now we consider the term mp— — q^KpA in (2) where mpHi 

oL • ' ' 



J 



Pj.Pp mjmp 



_ ms \Qj-Qp\ 

in (11). By our choice of ep, the closest distance between Jupiter and the planet 

is 0{m^l^). The interaction force between the two bodies is ^^^^^ — p{Qj — ^ 

\Qj-QpV 

1 /3 

whose modulus is mpO{mj ). Moreover, the force lasts only for a time interval of 

1/3 

order 0{m J ). Most of the time, the interaction force is 0{mjmp). This means 
the motion of the two bodies is only a small perturbation to the double two-body 
problem. Since pa^Qa are bounded when we average 5f along 71,72. Then the 

contribution of the term mp q\KpA to the averages 5fi is 0{mp'm j ). 

oLj 



mp 1 1 

Finally, we study the last term in (|2l) , say . We evaluate the average — / 

LI TAP Ti Jq Tap 

We know that 71 and 72 have almost elliptic orbits from theorem [Sl We choose the 
parametrization 7^ : (Lj + cos Kjt, 6i sin Kjt), due to expression ^51) (Recall the 
first appeared in (fsl)). We obtain 



dt. 



1 /"^^ 1 



Ti Jo rAP 



1 /" ' 1 

dt= — dt 

Ti Jo \J [xp — Li — ai cos KitY + {yp — h sin Hit)"^ 



1 /■ ' 

= — / [{xp - Lif + yp + {ai cos Kitf + {h sin Kitf 
Jo 

— 2ai{xp — Li) cos Kit — 2biyp sin Kit]~~^^'^dt. 

Denote [xp — Li)'^ + yp by rpi-, which is the distance between the planet and the 
Li Lagrangian point. The integral is taken w.r.t. t, while the position of the planet 
{xp,yp) is independent of t. So {xp,yp) has nothing to do with the integral, rp^. 
is much larger than the other terms under the square root. So we can use the Taylor 
expansion. 

= TT^^ — — / [{ai COS Kitf + {bi sin Kit)'^ -2ai{xp - Li) COS Kit 

rpLi ^rj,j^,Ti Jo 

-2b.yp sin K4]dt + 0(1 Akj = - + O (-^] . 

1 + a 

To summarize all the above computations, we obtain that when — t- 1 ^ — 

Ll 

1/3 — TTlp 

where a = dist(L2,(l — mj,0)) = OlmJ ), then 6/2 ~ blows up while 

rpL2 
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^/i ~ remains bounded by 0{mpm j ). 

This implies that /i(j^) — /i(i^) is nonconstant for this choice of ep. Due to the 
properties of analytic functions, fi{i') — fii^) is nonconstant for all other choices of 
ep. So we set ep = C such that our 4BP is far from collision. Then the proof of the 
lemma is completed. □ 



Appendix A. Delaunay variables 



The purpose of the appendix is to give a brief introduction to the Delaunay coordi- 
nates used in the paper. The materials could be found in [8j. For two-body problem 
of the form 

we know it is integrable in the Liouville- Arnold sense when H < 0. So we have the 
action-angle variables (L, i, G, g) to write the Hamiltonian can be written as 

H{L,£,G,g) = -^, 
The Hamiltonian equations are 

L = G = g = 0, 



{L,i,G,g) G T*T^ 



L3 



If we define the quantities: 

E: energy, M: angular momentum, e: eccentricity, a: semimajor, h: semiminor, 
then we have the following relations which endow the Delauney coordinates the 
physical and geometrical meanings. 



mk ' 



LG_ 
mk ' 



E 



k_ 
'2a'' 



M = G, e 




Moreover, g is the argument of periapsis and i is called the mean anomaly, and 
can be related to the polar angle ip through 



tan 



1 + e u 



esmu 



We also have the Kepler's law 



a 



1 



which relates the semimajor a and the 



r2 (27r)2 

period T of the ellipse. For two-body problem. Consider a body with position 
(^1,^2) and momentum (^1,^2)- we have the following formulas 



Q2 



a(cos u 



smu, 



Pi 

P2 



^mka ^/2- 



smn 



1 



e cos u 



cos u 



1 



ecostt 



aCNXIN XUE DEPARTMENT OF MATHEMATICS, UNIVERSITY OF MARYLAND, COLLEGE PARK, MD, 20742 

where u and I are related by it — e sin u = L 

Convert everything except u into Delaunay, we have the following 



(10) 



qi = {L'^/mk) |^cosu 
q2 = (LG/mk) sinu. 



-\ 1- 



G2 
L2 



Pi 



P2 



mk 



1 - 



mk 



smn 

1 ^ ' 
1 — ^ cos u 

COS u 



L2 



G2 

cosu 



Here g docs not enter because the argument of perihelion is chosen to be zero, g 



will enter if we rotate the (91,^2) and {pi,P2) using the matrix 



cos 5 
sing 



— sing 
cos 5 



Appendix B. Three-body problem in Delaunay variables 



For general three-body problem, the Hamiltonian is 

2 2 

Pi \ " mimj 



H3{qo,qi,q2,Po,Pi,P2) = ^w^ - XI 



2mi ^ \qi — qA 

where pi G M2 is the momentum, qi G M.'^ is the position, > is the mass. The 
first sum is the kinetic energy and the second sum is the mutual action potential. 
The system has 6 degrees of freedom. But it is also translation invariant and rotation 
invariant. We can eliminate 3 degrees of freedom after reducing the two invariance. 



B.l. The first step of symplectic transformation. We first eliminate the trans- 
lation invariance i.e. the momentum conservation. 



dpi A dqi = d{po + '^pi) Adqo + '^dpi A d{qi - go) := dPo A dQo + ^ dPj A dQi 
1=0 i=l 1=1 1=1 

We know dP^ = due to momentum conservation. 

{P1+P2Y . Pi . P2 momi mim2 mom2 



2mo 

1 1 

— + — 

771-0 "T-l 



+ 



Pj 

2mi 



+ 



Pi 



2m2 

m^mi 



\Q 



\Q 







+ 







m\m2 
ki - 92 1 

mo 777,2 



p2 
^2 



777om2 



+ 



P1.P2 

7770 



77717772 
IQ1-Q2I 



The Hamiltonian of three-body problem is now split into 2 two-body problems with 
a perturbation. 
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B.2. The second step of symplectic transformation. We convert the system 
into Delaunay coordiantes. Our symplectic form becomes 

2 2 2 

' d{Gi+G2)/\dgi+dG2/\d{g2-gi)+^dLiAdei 



y^^dPjAdQi = y ] dLjAdii+dGiAdg, 

i=l 1=1 



i=l 



The Hamiltonian becomes 

(11) = 



ki k2 , P1.P2 



2Ll 



2Ll 



+ 



mim2 



mo 



where ki 



\Qi-Q2\' 

i = 1,2 and we also need to convert Pi , Qi in the perturbation 



mo + mi 

into Delaunay variables. We notice that in the Hamiltonian, only the relative angle 
92— 91 appears. The rotation invariance implies the total angular momentum G1+G2 
is constant. Therefore we eliminate the term d{Gi + G2) Adgi in the symplectic form 
to obtain a system of three degrees of freedom. Now replace mo = 1 — mj,mi = 
"iJ;"2-2 = mp for our S-J-P system. So in our case, the corresponding geometric 
quantities. 



mj(l - mjfaj = mj(l - mjf, uj = 1, 



where c • 



m 



7-2 



kp 
Wp 



(1 — mj)'^mp 
1 — mj + mp 



Op := cap, 




mj + mp 




If we define mpHi 
are: 



Pl.P2 

mo 



(12) 



kj ^ 

L% oLj 



mim2 

IQi -Q2 
dHi 



in (11), then the Hamiltonian equations 

dHi 



1 + mp 



dHi 



'J 



9 J = mp 



dHi 



dGj 

kp 

n^'^'dTp^ 



Gj 
dHi 



gp = mp 



dHi 



Gp 



—mp- 
Lp 

■ —mp 



dHi 



Lj 



-mp 



-mp 



dHi 
dgp 



dHi 
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